Electrostatically confined Quantum Rings in bilayer Graphene 
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We propose a new system where electron and hole states are electrostatically confined into a 
quantum ring in bilayer graphene. These structures can be created by tuning the gap of the graphene 
bilayer using nanostructured gates or by position-dependent doping. The energy levels have a 
magnetic field (Bo) dependence that is strikingly distinct from that of usual semiconductor quantum 
rings. In particular, the eigenvalues are not invariant under a Bo — > —Bo transformation and, for 
a fixed total angular momentum index m, their field dependence is not parabolic, but displays two 
minima separated by a saddle point. The spectra also display several anti-crossings, which arise due 
to the overlap of gate-confined and magnetically-confined states. 
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The electronic properties of low-dimensional systems 
have long been an area of intense research. In re- 
cent years, a particularly interesting two-dimensional 
(2D) system, namely single layers of crystal carbon 
(graphene), has been obtained experimentally. The strik- 
ing mechanical and electronic properties of graphene have 
quickly been recognized and promise to lead to new appli- 
cations in electronic devices and sensors. These devices 
will benefit from the large charge carrier mobilities and 
long mean free paths at room temperatur e) 1 ' 2 ' 3 ' 4 ' 5 . It has 
been recognized that two coupled layers of graphene (bi- 
layer graphene - BG) is a very different material from 
graphene and also from graphite. The carrier spectrum 
of electrons in ideal BG is gapless and approximately 
parabolic at low energies around two points in the Bril- 
louin zone (K and K')&2. In the presence of a perpen- 
dicular electric field, however, the spectrum is found to 
display a gap, which can be tuned by varying the bias, or 
by chemical doping of the surface 8 . This tunable gap can 
then be exploited for the development of BG devices. In 
particular, the possibility of tailoring the energy gap in 
BG has raised the prospect of the creation of electrostati- 
cally defined quantum dots in bilayer graphene by means 
of a position-dependent doping or through the deposition 
of split gates on the BG surface^. 

A very important class of quantum devices consists of 
quantum rings. They have generated a lot of interest, 
especially because they allow the observation of quan- 
tum phase coherence effects on carrier transport such 
as in the Aharonov-BohmiS and Aharonov-Casher— ef- 
fects. Semiconductor-based quantum rings have been 
obtained experimentally 12 and, recently, quantum rings 
have also been studied on single-layer graphene, both the- 
oretically and experimentall y 13 ' 14 . The latter have been 
produced by lithographic techniques in which graphene 
nanoribbons or ring structures are carved from an oth- 
erwise defect-free surface. Such techniques permit the 
production of devices in the nanometer scale, but have 
the disadvantage of creating defects at the edges of the 
graphene-based structure which may reduce the overall 
performance, as well as raise difficulties for the theoreti- 
cal analysis of the device. These quantum rings can also 



be described as graphene flakes with a central antidot, 
and recent calculations indicate that localized edge states 
strongly affect the spectrum of these systems^. Addi- 
tionally, the specific type of edge (zig-zag versus arm- 
chair) was also found to influence their electronic prop- 
erties. 

To overcome the problems related to the edges of the 
quantum ring (i.e. defects and type of edge) we propose 
here an electrostatically defined graphene-based quan- 
tum ring (GQR). In contrast with previous quantum 
rings on graphene, the system considered here is cre- 
ated by the use of electrostatic potentials which induce 
a position-dependent gap, such that the low-energy elec- 
tron and hole states can be confined in a ring-shaped 
region of an otherwise ideal BG sheet. This type of con- 
finement is not possible in single layer graphene, due to 
the Klein tunneling effect^. In our proposed structure, 
the effects of the edges are not relevant since the BG 
sheet is assumed to be defect-free and the confinement is 
brought about by an external bias. Moreover, in contrast 
with structures carved by means of lithography, the ring 
parameters can be tuned by external fields, a feature that 
can be particularly relevant for the design of field con- 
figurable devices. Here, we obtain the energy levels and 
wavefunctions of the confined electron and hole states by 
numerically solving the Dirac equation and results are 
presented as function of ring radius, width and magnetic 
field. 

Bilayer graphene consists of two weakly, van der 
Waals coupled honeycomb sheets of covalent-bond car- 
bon atoms in a Bernal AB stacking. The system can be 
described in terms of four sublattices, labelled A, B (up- 
per layer) and C, D (lower layer). The A and C sites are 
coupled via a nearest-neighbor interlayer hopping term 
t. We employ the continuum model based on the Dirac 
equation, which is known to provide a realistic descrip- 
tion of graphene-based structures with dimensions that 
are much larger than the lattice parameter and that has 
been shown to give good agreement with experimental 
data. 

The Hamiltonian, in the vicinity of the K valley and 
in the presence of a magnetic field normal to the plane 
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FIG. 1: Schematic depiction of the potential profile for a bi- 
layer graphene quantum ring 



of the layer is given by 8 
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where t « 400 meV is the interlayer coupling term, 
7r = vp[(p x + eA x ) +i{p y + eA y )\, p is the momentum op- 
erator, A is the vector potential, AU = U\ — U 2 is the dif- 
ference of potential between the layers, U = (U\ + Z7 2 )/2 
is the average potential, and vf ~ 1 X 10 6 m/s is the 
Fermi velocity. In this work we neglect the small Zee- 
man splitting of the energy levels. 

The eigenstates of Eq. (1) are four-component spinors 
* = [tp A , , faj , ^d] t , where fa (j = A,B,C,D) 
are the envelope functions associated with the proba- 
bility amplitudes at the respective sublattice sites of 
the upper and lower graphene sheets. For a circular- 
symmetric, position-dependent potential such as the one 
described in Ref. [9J, the spinor components are if) a — 
Mpy m9 , i>B = Mp)^™- 1 ' 9 , faj = <f>c(pV me , and 
fao = 4>D(p)e l ^ m+1 * >S , where m is the angular momentum 
label. Using the symmetric gauge A = (0, B p/2, 0), the 
radial dependence of the spinor components is described, 
in dimcnsionless units, by 
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p/R, a = e — uo, uq = (u\ + U2)/2, and 



S = Aw/2, with Au = u% — u 2 - The energy, the po- 
tentials and the interlayer coupling strength are written 
in dimensionless units as f = E/Eq, Ui t 2 = Ui^/Eo, 
t' = t/Eo and Eq = Hvf/R, where p is the radial variable 
and R is the ring radius. The dimcnsionless parameter 
(3 = (eB/2h)R 2 , can be expressed as $/$o, where $ is 
the magnetic flux threading the ring and $0 = h/eis the 
quantum of magnetic flux. 

In this paper we solve Eq. (2) by the finite elements 
method for the following profile (see Fig. 1) : 







< r' < 1 



1 _ «. < r > < 1 + hl. 
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- < r' < 00 
(3) 

where w' — W/R, with W being the width of the ring. 
This potential describes a bilayer graphene quantum ring 
(GQR) of radius R and width W, in which both electron 
and holes arc confined by a tunable potential barrier Ub- 
As shown in Ref. [9], the wavefunctions given by Eq. (2) 
are eigenstates of the operator 
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with eigenvalue m, where I is the 2x2 identity matrix, 
L z is the angular momentum operator and a z is one of 
the Pauli matrices. The above equation can be rewritten 
as J z = L z + hr z +hS z , where the second term in the right 
hand side is a layer index operator, which is associated 
with the behavior of the system under inversion, whereas 
S z carries the information on the pseudospin components 
in each layer. It should be emphasized that in this sys- 
tem, L z does not commute with the Hamiltonian and is 
no longer quantized. 

Figure 2 shows the lowest energy levels of a GQR as 
function of ring radius R (upper panel), and width W 
(lower panels), for m = — 1 (a,d), m = (b,e), and m = 1 
(c,f), for W = 20 nm (a-c), R = 50 nm (d-f), U b = 150 
meV and B = 0. The results show a weak dependence 
on the ring radius and an interesting electron-hole asym- 
metry, which is due to the breaking of the inversion sym- 
metry by the external bias. However, the electron and 
hole energy levels are invariant under the transformation 
E(m) — > —E(—m). In Fig. 3, the lowest energy levels of 
a QR are shown as function of ring radius in the presence 
of a non-zero magnetic field, Bq = 5 T, for —3 < m < 3, 
W = 20 nm and Ub = 150 meV. In contrast with the 
zero field case, the results show a strong dependence on 
the ring radius, with the appearance of several crossings, 
as shown more clearly in the inset. 

The energy levels of electrons and holes for the K valley 
in a GQR are presented in Fig. 4 as function of exter- 
nal magnetic field Bo, for — 10<m<— 1 (red curves), 
rn = (green curves) and 1 < m < 10 (blue curves) 
with R = 50 nm, W = 20 nm and U b = 150 meV. In 
contrast with the results obtained from the Schrodinger 
equation, the spectra are not invariant under the trans- 
formation Bq — ► —Bq. Moreover, the individual branches 
display two local minima. The figures also show a set of 
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FIG. 2: Lowest energy levels of a graphene bilayer quantum 
ring as function of ring radius R (Panels: a, b, c) and width of 
the ring W (d, e, /). Parameters given in the text. 




FIG. 3: Energy levels of a bilayer graphene quantum ring as 
function of ring radius R. Parameters given in the text. 



branches with energies that decay as the field increases. 
These states correspond to Landau levels of the biased 
bilayer, for which the carrier confinement is mainly due 
to the magnetic field. By making a small gap in the con- 
finement region of the ring it can be evidenced that those 
branches are only weakly affected by a variation of the 
gap. We found that the electron and hole energies are 
related by E e (m, Bq) = —Eh(—m, —Bo). The spectrum 
for the K' valley can be obtained by setting Bq — ► — Bq, 
which shows that, as in the case of the graphene bilayer 
quantum dot£, the valley degeneracy is lifted. The four 
lower panels show the probability density of the points 
on the spectrum which are labeled by (a), (&),... in the 
upper panel. The vertical dashed lines indicate the edges 
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FIG. 4: Electron and hole energy states of a graphene bilayer 
ring-shaped quantum dot as function of external magnetic 
field Bq. Parameters given in the text. 




- Numerical results 
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FIG. 5: (a) Enlargement of the electron spectrum of Fig Q. 
(b) Number of transitions between the levels (circles) and AB 
oscillations for a ring with radius R and surface S = tyR 2 , 
S — iv(R ± W/2) 2 . (c) Ground state expectation values of 
L z /h, Tz/h and S z /h as function of the magnetic field. 
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of the ring. As seen, for the (e) and (f) points the prob- 
ability density has peaks in the internal region of the 
ring structure (p < R— W/2), whereas for the remaining 
points the carriers are mainly confined inside the ring. 
The (e) and (f) states belong to energy states that in 
the limit Bq — ■> oo approach the zero Landau level are 
shifted from E = due to the presence of the ring po- 
tential. An enlargement of the region of the spectrum 
is shown in Fig. 5(a). The yellow dots correspond to 
the location of the transitions of the ground state be- 
tween the levels which is related to the number of AB 
oscillations n = &/§o. Figure 5(b) shows the number of 
transitions between the levels (circles) and the AB oscil- 
lations for a zero-width ring with radius R and surface 
S = nR 2 , S = tt(R ± W/2) 2 . Notice that the numerical 
results obtained for the finite width ring follow closely the 
AB oscillations as expected for an ideal ring with radius 
equal to the one of the inner ring. 

The ground state expectation values of the L z /H, t z /H 
and S z /h operators are shown as function of the magnetic 
field in Fig. 5(c). The figure shows that the values of 

< L z > /H (black curve), < t z > jh (upper inset) and 

< S z > /h (lower inset) are not quantized, but, their 
sum, as in Eq. (4) (blue curve) is given by the to quantum 
number. Very different from the Schrodinger case we find 
that for Bq ks the quantum number to and the average 

< L z > are non-zero. 

The spectrum in Fig. 4 exhibits several anticrossings 
between different branches. An enlargement of one par- 
ticular anticrossing (between the first and second level of 
to = — 1) is displayed in Fig. 6(a). The corresponding 
values of < L z >, the average radial position < pjR >, 
as well as the probability density for the points (1), (2) 
and (3) are shown in panels (b), (c) and (d), respectively. 
The figure indicates that at the anti-crossing the eigen- 
state evolves from a configuration in which the electron is 
mainly confined in the center of the ring by the magnetic 
field, to a state in which the confinement is mainly caused 
by the electrostatic potential Ub- At the anti-crossing (2), 
the probability density plot shows peaks both in the cen- 
tral part of the structure, as well as in the ring, which 
characterizes the overlap of the magnetically and electro- 
statically confined states. 

In summary, we obtained the eigenstates and en- 
ergy levels of electrostatically confined quantum rings in 
graphene bilayers. This new approach differs from previ- 
ous studies of graphenc-based quantum rings which were 
realized through a lithographic cutting of the sample. 
The advantage of the present approach is that the con- 
finement can be tuned by varying the profile and depth 
of the confining potential, preventing issues as edge dis- 
order and the specific edge type. Numerical results were 
presented for a step-like confining potential, which can be 
expected to provide a good approximation of actual po- 
tential barriers, because of the large screening in bilayer 
graphene. We predict that the energy levels display a 
dependence on the magnetic field that is quite distinct 
from that of a conventional 2D electron gas. Moreover, 
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FIG. 6: (a) Enlargement around an anti-crossing point in the 
energy spectrum, (b) Expectation value of angular momen- 
tum < L z > jh and (c) radial average position < p/R > 
around the anti-crossing point shown in (a), (d) Probability 
density for the points (1), (2) and (3) shown in (a). 



the energy levels are not invariant under a Bq — > —Bq 
transformation, in contrast with usual ring structures. 
This is a consequence of the fact that the ring structure 
is produced by a gate that introduces an electric field 
and thus a preferential direction. On the other hand, 
the spectrum is found to be invariant under the trans- 
formation Bq — » — Bq, together with Ub — > —Ub and 
to — > —to. Alternatively, Bq — > — Bq and E — > —E 
transforms electron states into hole states and vice- versa. 
The spectra also present several anti-crossings at low en- 
ergies, which arise due to the overlap of gate-confined 
and magnetically-confined states. For a fixed total angu- 
lar momentum index to, the E(Bq) curves are no longer 
parabola, but show two minima separated by a saddle 
point. The existence of Aharonov-Bohm oscillations for 
both electrons and holes are still linked with flux quanti- 
zation through the ring. Throughout the calculation, we 
considered single electron states. However, the inclusion 
of Coulomb interaction between carriers may introduce 
interesting modifications in the spectrum^. The sys- 
tem can be realized experimentally by a suitable choice 
of nanostructured doping levels or with the application 
of nanostructured gates. The spectra may be investi- 
gated by, e.g. cyclotron resonance methods^ and quan- 
tum transport measurements. 
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